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Application of mean-field theory yields as main prediction that for symmetric 
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logarithm of the system size depending on the exponent of the degree distribution. 
Extensive numerical simulations contradict these results and clearly show that the 
mean-field assumption is not appropriate to describe this problem. 
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1. Introduction 

The zero temperature Glauber dynamics for tlie Ising model is one of the simplest 
ordering dynamics that can be devised for a set of spin-hke variables. Its sole physical 
ingredient is the tendency of each spin to align with the majority of its neighbors, and for 
this reason Glauber dynamics has been long investigated as a model for the ordering of 
magnets as well as of other systems where the evolution is driven by the minimization 
of interfacial energy without conservation constraints (0). Despite its simplicity, the 
zero temperature Glauber dynamics exhibits nontrivial phenomena already on regular 
lattices. In d = 1 its ordering process can be mapped to an ensemble of coalescing 
random walkers and, as a consequence, the mean ordering time grows with the system 
size N as A^^ (0). The same diffusive behavior occurs at T > in higher dimensions, 
but for strictly vanishing temperature there is a probability that the dynamics remains 
trapped in states not fully ordered (H; 0)- In ci = 2 the lack of global ordering occurs 
about 3 times out of 10, in the limit of large systems, and it is associated to the formation 
of two or more ordered stripes spanning the system along one direction. In three (and 
presumably higher) dimensions the probability of reaching the fully ordered state goes 
to zero as the system size increases: domains of opposite magnetization coexist, with 
some "blinkers" active at the boundaries 



The simplicity of its local interactions makes Glauber dynamics a natural model for 
the evolution of some social systems, whose dynamics can be defined in terms of agents 
characterized by a binary variable (for instance their opinion with respect to some issue) 
that is updated in response to the pressure of their peers. While for physical applications 
it is appropriate to study spins on regular dimensional lattices, the modelization of 
social systems immediately raises the question of the behavior of Glauber dynamics when 
the interaction pattern is disordered, i.e., agents live on the vertices of a network (j^. In 
this scenario, the nontrivial topological attributes of networks, such as the small-world 
property (0) and the scale- free property (3), are expected to play some relevant role. 
Steps in the direction of investigating these effects have been recently taken for other 
types of ordering dynamics, including the voter model llll l and the Axelrod 

model (Ui 



Among the first investigations of Glauber dynamics on networks we can mention 



the work of Boyer and Miramontes (|13[). who found that the long-range connections of a 
Watts-Strogatz network (j^ induce a pinning of the ordering process when domains reach 
a characteristic size, depending on the density of such connections. The Watts-Strogatz 
network is a rather special type of network, where, despite the small-world effect, the 
concept of a local ci-dimensional regular neighborhood of a site maintains a significance, 
at least for a small density of long-range connections. Other kinds of networks, such 
as the random (Erdos-Renyi) network (|1J) are instead not a perturbation of a regular 
lattice, and also in this case it has been shown that the zero temperature Glauber 
dynamics does not lead to full ordering (QEi). However, for networks of finite size 
A^, it occurs that some realizations of the process actually lead to complete order, while 
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others get trapped in an infinitely long-lived (meta)stable state. The probability that 
a run ends up in a disordered state turns out to grow with A^: order is not reached for 
asymptotically large (fl^. 

Very recently, Zhou and Lipowsky (jl7|) have presented a mean-field (MF) treatment 
of the zero temperature Glauber dynamics on generic uncorrelated complex networks 
whose degree distribution P{k), defined as the probability that a vertex is connected 
to other k vertices, decreases as a power-law, P{k) ~ k~'^, with a characteristic degree 
exponent 7. Within this approach, they find that a dynamic transition takes place for 
the particular value 7c = 5/2. For 7 > 5/2, the system reaches the ordered state in a 
time that diverges logarithmically with the system size. For 7 < 5/2, on the other hand, 
the ordering time is finite and independent from N. These results are claimed to be 
confirmed by numerical simulations. The results reported in Ref. (jl?! ) are remarkable 
in several respects. In the first place, they seem to contradict the difficulty observed to 
reach complete order generally observed for Glauber dynamics on disordered interaction 
patterns. Secondly, the transition in the dynamical properties occurs for a value of 
the exponent 7 that has not been associated to special behavior for other models on 
networks with a power-law degree distribution, in which transitions typically occur for 



7 = 3 ()18t Il9l : |2O0, below which the second moment of the degree distribution diverges 
and the network becomes scale- free (0). 

In this paper we present a refined study of the MF theory for the zero temperature 
Glauber dynamics presented in Ref. (0), recasting it more conveniently in a rate 
equation in the continuous time limit. In our analysis we take into account the finite size 
effects of the network in a more precise way, elucidating what happens at 7 = 5/2. In 
particular, we present a discussion of the differences shown by the dynamics depending 
on the initial conditions, that is, depending on whether the initial state is completely 
random and uncorrelated with zero average magnetization (symmetric initial state), 
or it corresponds to a state with a given initial, small, magnetization (asymmetric 
initial state). In the former case, we show that a proper description of the ordering 
process needs to take into account the initial diffusive contribution of the dynamics, not 
explicitly present in the MF equations for infinite network size. Our results confirm the 
MF prediction of a different behavior of the mean ordering time depending on whether 
7 is larger or smaller that 5/2. In order to check the predictions of our MF analysis, 
we report the results of extensive numerical simulations. These simulations prove that, 
for symmetric initial conditions, the dynamics behaves in a way incompatible with the 
MF predictions (as already observed numerically in other nonequilibrium dynamics on 
complex networks (j2ll )). showing no qualitative differences for 7 above or below 5/2, a 
value that seems to emerge as an artifact of an invalid MF theory. Our findings reconcile 
the results for Glauber dynamics on scale-free networks with those obtained on other 
types of structures and calls for approaches, beyond mean-field, able to capture the 
observed phenomenology. 
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2. Mean-field theory 

Let us consider an Ising model with spin variables a = ±1 located at the vertices of 
a network, which is fully characterized by the adjacency matrix Aij, taking the value 
Aij = 1 if the vertices i and j are connected by an edge, while Aij = otherwise. The 
local field hi(t) experienced by vertex i, and due to the spins located at its nearest 
neighbors, is thus given by 

h.it)=Y.A,a,{t). (1) 

j 

The Glauber dynamics at zero temperature is defined as follows: At each time step 
t we pick at random one of the vertices and we update the value of its spin according 
to the local field, that is, 

= + 1 if hi{t) > (2) 
ai= - I if hi{t) < (3) 
ai= ±1 with probability 1/2 if hi{t) = 0. (4) 

After each spin update, time is increased as t ^ t + 1/N, where N is the number of 
vertices in the network. 

Starting from a disordered uncorrelated configuration, in which each spin is taken to 
be positive at random with probability p, Glauber dynamics tends to make neighboring 
spins equal in sign. Thus, the main questions raised by this dynamics are whether all 
spins become equal in the long run and, if this is the case, how much time this total 
ordering requires, and how it depends on the initial fraction p of positive spins and on 
the system size A^. Particularly interesting in this respect is the symmetric initial state 
in which p = 1/2, since in this case the dynamics must first break the symmetry of the 
initial state in order to reach full order. 

To construct a MF theory for the zero temperature Glauber dynamics, we consider. 



following Ref . (|l7l ) , the average dynamical variable (t) , defined as the probability that 
a vertex of degree A; is in a +1 state. The evolution of the variable qk{t) depends on 
the average local field felt by the vertices of degree k, h^. In terms of this variable, and 
according to the rules defining the Glauber dynamics, we can write a rate equation for 
qk{t) in the continuous time approximation (which is valid for infinitely large network 
sizes) that takes the form 

= - qk{t)PToh[hk < 0] - qk{t)PToh[hk = 0]^ 

+ (1 - qk{t))Fioh[hk > 0] + (1 - qk{t))FToh[hk = 0]^, (5) 

where PToh[hk > 0], Prob[/ifc < 0], and PToh[hk = 0] are the normalized probabilities 
that the local field at the vertices of degree k is positive, negative or zero, respectively. 
Rearranging the terms in Eq. (0) we can write 

= -g,(t) + Prob[/^fc > 0] + ^Prob[/i, = 0]. (6) 
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In order to estimate the probabilities of the local field hk, the relevant quantity to 
consider is the probability Qk that an edge departing from a vertex of degree k points to 
a +1 spin. For a generic network, statistically characterized by its degree distribution 
P{k) and its degree correlations given by the conditional probability P{k'\k) that a 
vertex of degree k is connected to a vertex of degree k' fl22n . the probability that an 
edge from a vertex of degree k points to a +1 spin is given by 

Qu = Y.Pik'\k)qk, (7) 

k' 

that is, it is proportional to the probability that a vertex k is connected to a vertex of 
degree k' times the probability that this vertex is in a +1 state, averaged over all possible 
values of the degree k' . For random uncorrelated networks, the conditional probability 
takes the simplified form P{k'\k) = kP(k)/{k), and therefore the probability Qk is 
independent of k and can be written as (|l7l ) 
1 



(8) 



Assuming that the probability of having a positive spin at the end of an edge departing 
from a k vertex is independent of the values of the spins at the extremes of the other 
k — 1 edges (which corresponds to a MF assumption), we have that the probability that 
i edges from a k vertex point to positive spins is given by a binomial distribution. The 
local field at a A; vertex can be zero only if k is even and exactly half of its edges point 
to +1 spins. Thus 

' k 



Prob[/ifc = 0] 



k/2 



g'/'(i - Q) 



k/2 



(9) 



for k even. On the other hand, hk is positive when more than half of its edges point to 
+ 1 spins. Therefore, the probability of observing a positive local field is 

' ( k^ 



Prob[/ife > 0] 



k-l 



E 

l=\{k+l)/2\ \ ^ / 

where [x] is the smallest integer larger than or equal to x 
the rate equation for the probabilities as 

dqk{t) 



dt 



-qk{t) + ^k{Q), 



(10) 

In this way, we can write 

(11) 



e,k/2 



k 



where fllTI ) 

k 

MQ)= E 

l=\k/2-] 

By introducing Eq. ()11|) into the definition Eq. 
the quantity Q, namely 

g = -g + v]>(g), 

with 

\'^/ k 



Q'il - Q)''-'. (12) 
3) we obtain a closed rate equation for 

(13) 
(14) 
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It is useful to write the dynamics in terms of the new variable y = Q — 1/2. In the 
initial state, where all spins point up or down with probability p and 1 — p, respectively, 
we have Q{t = 0) = p and therefore y{t = 0) = p — 1/2. In the symmetric disordered 
state, with p = 1/2, we have y{t = 0) = 0, while the ordered states with all spins up 
or down correspond to ?/ = 1/2 and y = —1/2, respectively. Thus, y plays the role of a 
convenient order parameter. The rate equation for the quantity y can be written as 

y = -y -^ + i^{y), (is) 

where we have defined the new function 

*fe)-*(i+.)=i:^*'=G+.). (16) 

The function ip{y) is equal to 1/2 for y = and goes symmetrically to 1 [0] for 
y = 1/2 [y = —1/2]. In terms of this variable, we can study the time that it takes for 
the system to become ordered by defining the ordering time tord such that \y(tord) \ = yp, 
where yp is an arbitrary chosen positive value such that \p — l/2\ < yp < 1/2. 

^From the explicit evaluation of ip{y), presented in the Appendix, it turns out that, 
for an uncorrelated network with degree distribution P{k) ~ k~'^, there are different 
limiting regimes for different values of y and the degree exponent 7. For simplicity we 
consider only the case y > 0; analogous considerations apply in the symmetric case 
y<0. 

2 < 7 < 5/2: The behavior of ip{y) can be approximated, depending on the value 
of y and of the minimum degree ko and degree cut-off kc of the network in the following 
form: 

ijiy) - 1/2 ~ C2i^)k2-'kl/'-yy, y<y, 

i^{y)~i/2^c,{^)krV^''''\ yc<y<yNL (17) 

ij{y) - 1/2 ~ 1/2, y > y^L 

where y^ = k'^^^'^ and y^L = ^o^''^- 

Inserting these forms into the dynamic equation Eq. (fTH|l . we obtain the temporal 
evolution of the order parameter y. For y < y^, there is a regime with y growing 
exponentially with time, i.e. 

y{t) = yrexp[{c2{j)kr'kl/'-^ - l){t - t^)], (18) 

with an exponential factor that depends on the network cut-off kc. Here we have denoted 
as tr and yp the values of t and y at the beginning of this regime. For reasons that will 
be clear below, these may or may not coincide with the initial values t = and y{0). 
For yc < y < yNL, on the other hand, the order parameter grows instead as a power-law, 

yit) = [yl-'^'-'' + c,{i)kr\t - g] , (19) 

where tc is the time at which the order parameter reaches the value yc. Finally, for 
y > yNL, i^iy) is approximately equal to 1 and the dynamic equation is = 1/2 — so 
that 

y{t) = l+{yNL-l)e~^'''--\ (20) 
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where t^qi is the time at which the value is reached, that is, y{tNi) = Unl- 

7 > 5/2: In this case Uc does not play any role; the function il){y) is independent 
of kc and it can be approximated for y < i/atl as 

V'(y)~l/2 + ci(7)A;yV (21) 

The behavior y is thus again exponential, but now with a prefactor independent of k^ 
namely, 

y{t) = yTeM{ci{l)ky^ - - ^t)], (22) 

For y > yNL, ip{y) the scenario is the same of the previous case, with an exponential 
approach to the asymptotic value y = 1/2. 

As a consequence of the previous expressions, the total ordering time has, in the 
most general case, four contributions, namely, 

tord = tx + Ate + ^^NL + ^tord-, (23) 

where At^ = tc-tT, AtNL = ^nl - tc, and Atord = tord - ^nl- 

The simplest term is the last one, Atord, that, from Eq. (jSHI), is given by 

1 - 2yNL 
_l-2yF _ 

It is independent of A^, p, and, in practice, also of k^. The other terms depend instead 
on the value of 7. 

For 7 > 5/2, Ate and At^vL are merged in a single term, which is, from Eq. 

At, + AtNL^ko^'^\^ ^ 
For 2 < 7 < 5/2, instead, we obtain from Eq. (jl9|) 

1-2(7-2) _ 1-2(7-2) 

At^L - , (26) 



Atord ~ In 



(24) 



(25) 



k^ 



while from Eq. (fTHj) we have 



At, ^ kl-'<kl-'''^ In 



Vc 



(27) 

Let us now discuss the additional term tr- The deterministic equation of 
motion ()15p holds strictly only in the limit of infinite network size. For finite an 
additional random term appears, giving a diffusive contribution to the dynamics of y 
that competes with the deterministic drift described by Eq. (fT3j) . As it will be shown 
below, this contribution vanishes as A^ diverges, so that it can be generally neglected 
with respect to the drift, unless y{0) = 0, i.e., for symmetric initial conditions. In such 
a case also the drift vanishes, and the point ?/ = is an unstable equilibrium point, that 
can be characterized by a potential V{y) = — P{ip{y) — y — ^/2)dy with a parabolic 
shape for very small y. Under the deterministic MF dynamics a system placed at ?/ = 
will stay there forever, and the diffusive term is thus needed to provide the perturbation 



Zero temperature Glauber dynamics on complex networks 



8 



necessary to take the system out of the initial equihbrium point. The interplay between 
drift and diffusion determines the time tx at which drift starts to dominate. 

The origin of this diffusive term is easy to understand: As discussed at the beginning 
of this Section, the Glauber dynamics at zero temperature is defined by means of a 
sequential updating, in which at each time interval At = a single spin is flipped, 
causing an average increment of the dynamic variable AQ = 1/N. The initial condition 
y = 0, Q = 1/2, corresponds to a state in which every site is randomly assigned a +1 
or a —1 spin with probability p = 1/2. Therefore, the initial evolution is dominated 
by a random flipping of spins, driven by the value of the local field hi, that induces 
a diffusive motion of the variable y until a nonzero value is reached, large enough to 
drive the dynamics out of the unstable equilibrium point. To be more concrete, when 
?/ ~ 0, or Q — 1/2, the drift term appearing in the equation of motion is very small. In 
particular, the deterministic drift velocity is given by 

v = {l- Q)^{Q) - Q[l - = ij{y) - 1/2 - y. (28) 

The diffusive component, on the other hand, has a diffusive constant that can be 
estimated as 

D = (Ag)V(2At) = 1/{2N). (29) 

For small values of y, the drift grows with y, while diffusion is independent of it. Starting 
with y{0) ~ 0, for small \y\ diffusion prevails and the drift can be neglected. The diffusion 
leads to increasing fluctuations, and after some time a fluctuation leads to a threshold 
value sufficiently large for the bias to start prevailing, and drive the variable far 
from the equilibrium point. From this time on, bias dominates and the diffusive term 
in the equation of motion becomes negligible. 

The value of yT can be determined as follows: When the system is in the position 
y, the time to needed to go back diffusively to is given by \y\ = y/Dtn- During this 
time interval the bias will induce a drift \v{y)\t£, = \v{y)\y'^/ D. If the drift is smaller 
than \y\ the system is able to go back to 0, and diffusion dominates. If the drift is larger 
than \y\ then bias dominates. The condition that defines yr is then 

vivr)^ = IvtI- (30) 
Inserting Eq. (PT|) into Eq. (jHUj) we obtain, for 7 > 5/2, 

lyrl ^ N~'/%'/\ (31) 
while the time spent diffusing is 

tT - f - k~'^'. (32) 

Notice that this time does not depend on A^. On the other hand, for 2 < 7 < 5/2, we 
obtain, using the scaling of the cut-off with the network size, given by k^. ~ N^^"^ for 
uncorrelated scale-free networks (jisl). 
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Figure 1. Sketch of the behavior of the evolution of y{t) as a function of time for 
2 < 7 < 5/2. Both axes are on logarithmic scale. 



and 



a time that goes to zero when increasing A^. We are now in the position to summarize 
the MF estimates for the ordering time. 

For asymmetric initial conditions (p 7^ 1/2, ?/(0) > 0) we can set = and 
Ut = 1/(0). In this case, we find at leading order for 7 > 5/2 

tord ^ const + kQ^^"^ \ny{0), (35) 

while for 2<7<5/2we obtain 

" N 



tord ~ const + A;'"^Ar(^-5/2)/2 



y(o) 



(36) 



As we can observe, for finite y{0) the ordering time is finite in the infinite network size 
limit. For finite networks, it exhibits a decreasing size correction in the case 2 < 7 < 5/2. 
For symmetric initial conditions, y{0) = 0, and 7 > 5/2 we have, at leading order, 

tord ~ const + ko^^^ In N, (37) 

that is, the ordering time diverges logarithmically with the network size A^. For 
2<7<5/2, on the other hand, we obtain 

tord ~ const + A;^-^ Ar(^-5/2)/2 ^_ ^gg^. 

The ordering time is now a decreasing function of A^ that tends to a constant value in 
the limit A^ — 00. For illustration purposes a sketch of the different temporal scales 
and the related behaviors of y is depicted in Fig. [H 



3. Numerical results 



In this Section we compare the theoretical results derived above within the MF formalism 
with extensive numerical simulations of the Glauber dynamics at zero temperature. 
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Figure 2. Plot of the probability that a run starting from the symmetric initial state 
?/(0) = does not lead to complete order as a function of TV, for networks with various 
values of 7 and k^. 



performed on uncorrelated networks with a power-law degree distribution. The networks 
used have been generated using the uncorrelated configuration model (UCM) defined in 
Ref. (j2J). The network is built as follows: We start by assigning to each vertex z in a set 
of initially disconnected vertices a random degree fcj, extracted from the probability 
distribution P{k), subject to the constraint kg < ki < k^ = N^/"^ and X^i even. 
Afterward, the network is constructed by randomly connecting the vertices with 
edges, respecting the preassigned degrees and avoiding multiple and self-connections 
between vertices. Using this algorithm, it is possible to create random networks with any 
preassigned degree distribution. For the case of networks with a degree distribution given 
by a power law, P{k) ~ k~^, the restriction k < iV^/^ for any degree exponent guarantees 
that the networks generated in this way are completely uncorrelated (ji^. Notice that 
if such a restriction is not enforced, correlations unavoidably arise for 7 < 3 (ji^. 

We start considering the case 7 = 4 as a representative of networks with 7 > 3. In 
this case the connectivity pattern does not show strong degree fiuctuations ((A;^) < 00 
for — > 00) and one does not expect to observe the special effects induced by the 
presence of hubs, relevant for other dynamical processes on networks. MF theory 
predicts the ordering time to be given by Eq. (jH^ for y{0) > and by Eq. (jHTj) for 
y{0) = 0. However, even in this simple case, we find that the observed phenomenology 
is in striking disagreement with the expected theoretical results. 

The first, fundamental, discrepancy lies in the very existence of order in the system 
for large times. At the MF level, all instances of the dynamics order sooner or later. From 
the numerical evidence, however, we conclude that networks starting from a symmetric 
initial state do not get ordered in the large A^ limit. To reach this conclusion, we 
have performed a large number of runs for several values of A^, starting from y{0) = 0. 
We observe that some of them do not lead to a fully ordered state, remaining instead 
trapped in some (meta)stable states that last forever. Fig. |2] reports the fraction of 
these runs that do not lead to complete ordering. After a minimum for a network size 
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Figure 3. Average value of the ordering time tord: for 7 = 4 as a function of the 
network size, for fixed y(0) = 0.01 and different values of fco. Inset: Asymptotic 
ordering time as a function of the initial y(0). 



that apparently depends on the value of /cq, this fraction grows with N . Extrapolating 
for large A^, we conclude that an infinite network never gets ordered. This absence 
of ordering of the Glauber dynamics a T = on complex networks had already been 
pointed out in Ref. JigI ). Fig. |21 demonstrates also that the lack of ordering occurs for all 
values of 7, either larger or smaller that 5/2. From this figure we can also observe that 
the fraction of nonordering runs decreases with increasing fco? in such a way that for 
very large values of fco and moderate values of N , this fraction is apparently zero (i.e. all 
realizations get ordered). This fact, however, is a trivial consequence of the nonordering 
probability becoming too small for the number of runs that can be performed. In 
Ref. ()12i), where large values of fco were used, this absence of ordering was not noticed. 

In view of this result, the very notion of an ordering time becomes somewhat fuzzy, 
in the sense that an unrestricted average over all runs will yield an infinite ordering 
time. For practical purposes, nevertheless, we can define an operational ordering time 
by performing an average over the restricted ensemble of those runs that actually order. 
Thus we have computed the average ordering time, defined as the average time tord 
needed for the ordering parameter y to reach the arbitrary value \yp\ = 1/4 (|l7|). In 
Fig. IHl we report the values of tord obtained for networks with 7 = 4, starting from 
nonsymmetric initial conditions y{0) > and different values of fco, as a function of the 
network size A^. As predicted from the MF result Eq. (jH^j) . we observe that tord becomes 
independent of A^, for sufficiently large network sizes. In the inset of Fig. IHl we plot the 
asymptotic value of tord for the largest network size considered A^ = 10^, as a function 
of y{0). Here we can observe that the MF prediction tord ~ — ln?/(0), Eq. (jHSI), seems to 
be fulfilled, at least for sufficiently large values of fco and y{0). For small values of y{0) 
and fco the convergence of tord with A^ is rather slow, a fact that hides the asymptotic 
logarithmic behavior. 

In Fig. |3]we show the results obtained for the ordering time in networks with 7 = 4 
and symmetric initial conditions y{0) = 0, for different values of the minimum degree 
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kf). As we can see from this plot, the average ordering time tord seems to grow with 
as a power- law, with effective exponents depending on the particular value of k^. For 
large values of /cq, the exponent is so small that the behavior could be confused with a 
logarithmic one. However, for small /cq, the power-law growth is quite evident, in strong 
disagreement with the MF prediction, Eq. (jHTf) . 

We have performed this same kind of analysis for two other values of the degree 
exponent, namely 7 = 2.75 and 7 = 2.25, in order to check the possible effects of a 
strongly inhomogeneous degree distribution, with diverging second moment, and also 
to asses whether the values of 7 smaller than 5/2 affect the behavior of the ordering 
time, as predicted by MF theory. Results for asymmetric initial conditions y(0) > 
for 7 = 2.25 are reported in Fig. El In this case, we find that the convergence to the 
asymptotic value of the ordering time for large N is extremely slow. As a consequence, 
it is not possible to make a definite statement about the validity of MF based on the 
present numerical evidence. For the largest A^ that could be reached, the behavior of tord 
for small ?/(0) does not seem to be logarithmic, inset in Fig. |S1 in apparent contradiction 
with Eq. (jHEI)- However, it is possible that this discrepancy could disappear if larger 
values of A" could be considered. 

Much clearer is the situation for the symmetric initial condition |/(0) = 0: the 
behavior is qualitatively the same (and different from MF) independently of the value 
of 7. We have already seen in Fig. |2l that the absence of ordering for large A^ holds for 
any value of 7. But the similarity is stronger when the ordering time tord is measured as 
a function of A^. As displayed in Figs. El and [71 tord diverges at large A^ as a power-law, 
as opposed to the logarithmic growth and the constant behavior predicted by MF for 
7 = 2.75 and 7 = 2.25, respectively. 

The conclusion of this numerical analysis is that MF theory is unable to account 
for the behavior of the ordering time of the Glauber dynamics at zero temperature with 
symmetric initial conditions. The physical reason for this failure can be understood by 
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Figure 5. Average value of the ordering time iord, for 7 = 2.25 as a function of 
the network size, for fixed ?/(0) = 0.01 and different values of fco. Inset: Asymptotic 
ordering time as a function of the initial y(0). 
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Figure 6. Average value of the ordering time tord, for 7 = 2.75, as a function of the 
network size, for symmetric initial conditions 2/(0) = and two values of the minimum 
degree Uq. 



considering in more in detail the dynamical evolution of the system. 

The symmetric initial state is created by randomly placing +1 and —1 spins on the 
network with probability p = 1/2. This procedure is akin to a percolation process (ji^ of 
both the positive and negative spins with the same probability p. In random networks, 
the percolation threshold is given by (|2Cl ) 

1 , , 

= MTT' 

In scale-free networks with 7 < 3, the degree fluctuations diverge, (fc^) — > 00, and Pc 
in the limit of large network sizes. In the case of homogeneous networks, with finite 
(fc^), we observe that pc ~ l/(^); and thus it becomes small for large average degree 
or fco- Hence for reasonable values of ko, all networks exhibit a very small percolation 
threshold. Therefore, symmetric initial conditions give rise to two giant components. 
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Figure 7. Average value of the ordering time to^d, for 7 = 2.25, as a function of 
the network size, for symmetric initial conditions 2/(0) ~ and three values of the 
minimum degree feg- 




Figure 8. Temporal evolution of the size of the giant components {G+ and G-), 
the size of the boundary B and the order parameter y for a run on a network with 
N — 10"', 7 = 4, fcg = 5 with symmetric initial conditions. 

G-i- and G-, with approximate size N/2, each one corresponding to a connected set of 
vertices occupied by +1 or —1 spins, respectively. In Fig. IHlwe plot, as a function of 
time and in a single run, the relative size of the giant components of +1 spins, Gj^{t), 
and —1 spins, G-{t), and the relative size of the boundary B{t), defined as the fraction 
of edges with spins of different sign at the ends. For illustration, we also plot the value 
of the order parameter y{t) corresponding to that particular ordering run. From this 
Figure we can observe that the dynamics is governed by the competition of the giant 
components of different signs, the (?+ growing at the expense of the (?_. It is noteworthy 
that during most of the evolution the boundary occupies a constant and sizeable fraction 
of all the edges. In the case shown in the Figure the competition ends abruptly with 
the positive giant component invading the whole system. In some cases instead, 
the competition ends up in a stalemate: the system reaches a stationary configuration 
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Figure 9. Plot of the probability that the local field of a vertex with fco = 5 links 
is positive. Data are for a system with 7 = 4 and N = 100000. The meaning of the 
different curves is explained in the text. 



with the two coexisting giant components basically frozen with only some spins of even 
degree freely flipping back and forth at the boundary. 

In the initial state B{t = 0) = 2p{l —p). After a short transient B{t) reaches the 
plateau. During this short time interval strong correlations build up in the system, in 
the sense that +1 spins become surrounded with high probability with a majority of 
positive spins, while most neighbors of a —1 spin are almost certainly negative. This fact 
can be quantitatively measured by computing the mean probability that a vertex with a 
given spin a has a positive local field h. This kind of measurement is reported in Fig. IHl 
for homogeneous networks with degree exponent 7 = 4, where we display the probability 
that the local field of a site of degree ko is positive. Within the mean-field treatment, 
this quantity is given by the function $^(,(1/2 + y), Eq. (O, plotted as a dashed line. 
The solid line in the middle is the numerical evaluation of the same quantity. The other 
two lines are instead numerical evaluations of the same probability, but averaged only 
over all sites with Sj = 1 (top) and Sj = — 1 (bottom). The middle solid line (that is 
the average of the top and bottom lines) is quite similar to the analytical result, but 
the other two are completely different. The mean-field treatment assumes that the local 
field is completely uncorrelated from the value of the spin, and depends at most on the 
degree of the vertex in which the spin is located. In the real system, however, during 
virtually all the dynamics, a positive spin has almost surely a positive local field, while 
a negative spin has a negative field. These strong correlations, completely missed by the 
MF approach, are at the root of the nontrivial behavior shown by the zero temperature 
Glauber dynamics. 

4. Conclusions 

In this paper we have investigated the ordering process of the zero temperature Glauber 
dynamics when Ising spins are placed on the vertices of an uncorrelated network with 
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degree distribution decreasing as a power-law of exponent 7. We have first developed 
a mean-field theor y fo r the problem, expanding and improving on the earlier approach 
presented in Ref. (llTl ). Within the MF theory we identify a suitable order parameter 
and compute in detail the time needed for the system to order. While for asymmetric 
initial conditions the ordering time is finite, for symmetric initial conditions there is a 
transition: tord is finite for 7 < 5/2, while it diverges logarithmically with the system 
size for 7 > 5/2. 

The validity of these analytical results has been checked by means of extensive 
numerical simulations. While for asymmetric initial conditions we find an apparently 
reasonable agreement with MF theory, when the probability of a positive initial spin is 
p = 1/2, the behavior is at odds with the MF results and highly nontrivial. In particular, 
it turns out that there is a finite probability for a run to get trapped forever in a 
stationary disordered state. This probability grows with so that complete ordering is 
not reached for asymptotically large A^. This conclusion, in striking disagreement with 
the MF predictions, holds for all values of 7. If we compute the ordering time restricted 
only to those runs that actually lead to full order, other unexpected results arise. For 
large and independently of the degree exponent, tord appears to diverge as a power- 
law of the system size, the effective exponent being a function of 7 and the smallest 
degree k^. One of the main conclusions of our work is, therefore, that mean-field theory 
does not provide a correct description of the ordering dynamics of the Glauber model 
at zero temperature. In this sense, the peculiar exponent 7c = 5/2 turns out to be an 
artifact produced by a MF theory that is not capable to describe the system, and it 
does not bear apparently any relevance for the true behavior of the Glauber dynamics. 

The failure of the MF theory for symmetric initial conditions is reminiscent of 
similar observations made on other dynamics on scale-free networks 1)2 ih . In the present 
case, this failure can be traced back to the breakdown of the assumption that the local 
field of a spin is independent of its value. This breakdown is explicitly shown in the 
formation of two giant components of positive and negative spins, competing through 
an extensive boundary, which induces the presence of strong correlations. 

It is important to remark that the nontrivial phenomenology uncovered in this 
paper is not a consequence of the scale-free property of the network, since it appears 
already for 7 = 4. The nontriviality of this problem is thus totally unrelated from 
the heterogeneous or homogeneous nature of the substrate. Somewhat surprisingly, the 
behavior is qualitatively the same for any 7 > 2: In this case a scale-free topology seems 
not to have any relevant effects. The key to a deeper understanding of this problem lies 
likely in a detailed analysis of the properties of the two giant clusters competing in the 
system and of the boundary separating them. Along with an investigation of the effect 
of temperature in the problem, this remains as an interesting open issue for future work. 
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Appendix A. 

In this Appendix we derive analytically the asymptotic behavior of the function 



. {k) ^2 . 

for finite networks, in which the the degree k is restricted to the range A; G [fco, ^c], 
being the smallest degree in the network and kc the degree cut-off or maximum degree 
present in the network. In the continuous degree approximation, the degree distribution 
takes in this case the form 

The function $^(1/2 + y) has the generic form, neglecting the Kronecker symbol, 

k 

^k{l/2 + y)= ^ PkAv). (A.3) 

l=k/2 

where 

Pm(i/)=(^)(1/2 + Z/)'(1/2-1/)'=-^ (A.4) 

For large values of k, the binomial distribution Eq. ()A.4|) can be approximated by a 
normal distribution, 

Pk,i{y) ^ ^= exp (-(/ - mfl2a^) , (A.5) 



with 



m = Y.lpuM = k{l/2 + y), (A.6) 
I 

= E ^'pM -^^ = Ml - %')/4. (A.7) 

If we insert this form into the definition of $^(1/2 + y\ and substitute the sum over I 
by an integral, we obtain 



where erf(x) is the error function (26). For constant y^k and large fc, the argument 
of the second term in Eq. ()A.8|) becomes large, and we can approximate this term by 
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unity. In the same limit, the y contribution in the denominator of the argument in the 
first term of this equation is neghgible, so that 

$fc(l/2 + y) = l + ^erf {yV2k) . (A.9) 

In Fig. lAll we check the vahdity of this analytical expression by plotting it agains the 
result of the numerical evaluation of Eq. ()A.3|) as a function of the rescaled variable 

Inserting the functional form Eq. ()A.9|) into the definition of we obtain, within 
the continuous degree approximation, and neglecting terms of order {kc/kQY~"' , 

ilj{y) ^ i + ~ f^^ fc-^+^erf (yV2k) dk. (A. 10) 

2 2, J ko 

The behavior of this function has different asymptotic regimes, depending on the value 
oiy. In the first place, if \/2k^y is small, then the argument in the error function is small 
for all values of k, and we can substitute it by its linear approximation, erf (x) ~ 2x/ y^. 



Table Al. Functional forms of ip{y) — 1/2 as a function of y for the different ranges 
of 7. 





2 < 7 < 5/2 


7 > 5/2 


Vc-^y vnl 
y > yNL 


c,{^)kr'kl/'-^y 

C3ii)krv^^-'^ 
1/2 


ci{l)kl^^y 

ci{i)ky^y 

1/2 
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Then, for y < = we have 



2 0r Jfco 

^.^^lil^.^^^^,;^. (A,12) 

2 Vtt 5/2 - 7 

As we can observe, the behavior of ipiv) for V < He depends now also on the value of the 
degree exponent 7. For 7 > 5/2, the contribution of vanishes in the limit of large kc, 

1 /2 

and the function il){y) takes the functional form tp{y) ~ 1/2 + 01(7)^0 U- For 7 < 5/2 
on the other hand, the most relevant contribution comes from the degree cut-off, and 
we have ip{y) ~ l/2 + C2(7)/co~^^c^^~'''?/- Here 01(7) and 02(7) are constants that depend 
only on the degree exponent 7. 

When y > y^, we must take into account the full functional form of the error 
function in order to evaluate Eq. (jA.lOj) . In this case we have that the argument of the 
error function is large at the upper limit of the integral. Given that the error function 
saturates to one for large argument, we can simplify the integral extending its upper 
limit to infinity. Therefore, we have that, in this range of y values, 

7/^(7/) ^ ^ + ~ k-^+hif (y\/2k) dk. (A. 13) 

This last integral can be expressed analytically in terms of the incomplete Gamma 



. . —1/2 

function ()26f l . Performing an expansion of the resulting expression for y <^ yNL = k^ 



we obtain again two limiting regimes, depending on the degree exponent. For 7 > 5/2, 

1 /2 

the leading term is again ip{y) ~ 1/2 + ci(7)A;o y, that is, we obtain the same linear 
behavior as in the case y < yc- For 2 < 7 < 5/2 instead, the integration over k brings 
a singular behavior, so that we obtain ip{y) ~ 1/2 + 03(7)^0"^?/^*^'^^^''. Here 03(7) is 
again a constant that depends only on the degree exponent 7. Finally, for y ^ yjqi the 
argument in the error function in Eq. ()A.13j) is large for all k, and we can approximate 
erf(x) ~ 1. In this case, V'(y) can be taken to be equal to 1 for any value of 7. The 
behavior of the function ip{y) is summarized in Table \K\\ 
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